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1. Abstract 

An algorithm is described that localizes a set of 
simultaneously energized coils using MEG 
detectors. The algorithm is based on a 
mathematical model describing the coils as 
stationary magnetic dipoles, of which the source 
time functions are known. It is assumed that the 
source time functions are orthogonal. It is shown 
how this a priori knowledge can be used to improve 
the efficiency of the algorithm. 

The method is useful for continuous or intermittent 
head position registration during long MEG 
sessions, to registrate MR and MEG data, to 
localize EEG electrodes attached to the head, when 
EEG and MEG are recorded simultaneously, and to 
digitize the head shape. Experimental data shows 
that the localization error can be smaller than 
2 mm. 

2, Introduction 

Compared to Electro Encephalography (EEG) 
Magneto Encepalography (MEG) has the practical 
advantage that no electrodes need to be attached to 
the subject’s head. In particular, when many 
cha nn els are recorded for source localization 
studies, this is an important practical advantage. 
The drawback is, however, that the subject can 
move his head with respect to the MEG sensors 
during the registration. For this reason it has been 
proposed [1,3] to attach coils on the head which are 
energized with a sine function. When these coils 
are localized by fitting a magnetic dipole function 
on the resulting MEG signals, the position of the 
head with respect to the sensor system can be 
determined. When this procedure is performed 
repeatedly, the movement of the head can be 
recorded and the data can be corrected using 
interpolation techniques. 

In this paper a new variant of this head localization 
technique is described. Instead of using coil sets of 
three orthogonal coils [1,3] simple uncalibrated 
coils are used and the source time function of the 
currents are not restricted to sine waves, but can 
have any shape. Furthermore, the coils can be 
energized simultaneously. The proposed technique 


of localizing energized coils also has other 
applications. It can be used to project current 
dipoles, localized with MEG, onto an MR scan [2], 
Furthermore, when three coils are fixed to the 
subject’s head, a fourth coil can be used to 
determine the shape of the head or the positions of 
EEG electrodes, when EEG and MEG are recorded 
simultaneously. 

3. Methods 

When the current activation function of coil k is 
denoted by s k ((j) (or s kj ), and the resulting magnetic 
induction recorded from sensor i is denoted as by , 
we have 

i if = £^PGi.n i ;m k ,y k ) ,• i= 0,...,/-1 and j=0,...,J—l. 

Here, K is the number of coils, I is the number of 
sensors, J is the number of time samples and (3(x„ 
n,; m kj y k ) is the magnetic induction caused by the 
k -th coil on sensor i. Furthermore, Xi and n, are the 
position and orientation of the /-th MEG sensor and 
m k and y k represents the position and orientation of 
the k -th coil. 

When it is assumed that the background noise is 
uncorrelated, the optimal estimator of the coils is 
obtained by minimizing H, the sum of squared 
differences between the model and the data. 
Herewith, it is assumed that the source time 
functions of the currents applied onto the coils are 
known, up to a constant factor, which appears as 
the magnitude of m k in the equations. It is found 
that 

H= min 
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with 
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Note that since the source time functions sy are 
known a priority, the cross correlations of the 



source time functions Skua and the cross 
correlations between the source time functions and 
the MEG signals c± can be computed in advance. 
As a result, no summations over time are required 
in the evaluation of the cost function and the speed 
of the algorithm is (almost) independent of J. 


The cost function is split in separate terms, which 
depend on one coil. Therefore, when the source 
time functions are mutually orthogonal, solving a 
minimization problem for each of the coils 
independently solves the minimization problem for 
all coils. 


3.1 The optimal moment parameters 


Adding a pointing tip to the fourth coil 


When minimizing the function H, one can use the 
fact that the coil moments are linear parameters. 
This implies that when the gradient of H is set to 
zero, a linear system of equations in m T = (m 0 T , 
mi T ,...,m K -i T ), is obtained: 
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or 


Am = b. 

With the same notation the non-linear part of the 
minimization can be expressed as: 

H= min - m T b + Ho • 

y, 

With this expression for H, additional summ a tions 
over i are avoided. 


3.2 Orthogonal time functions 

When the applied currents are devised in such a 
way that the source time functions are mutually 
orthogonal the matrix A becomes 
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and the vector b can be expressed as b T = (b 0 T ,..., 
bit-i T ). Under these conditions one obtains for the 
residual function H 
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One of the applications of our coil localization 
algorithm is to attach three coils to the head, e.g. in 
front of both ears and at the nasion, and to use the 
fourth coil to localize EEG electrodes or to 
measure the shape of the head. Since the dimension 
of the coil (ca. 5 mm) is larger than the required 
and expected accuracy of the localization algorithm 
it is useful to add a pointing tip onto the fourth coil, 
which points to the measurement point. When the 
length l of the tip is known and one has localized 
the center y and the orientation m = m/|m| of the 
coil, one can compute the position t of the tip using 
t=y+/m . 



Figure 1. When a pointing tip is attached to the 
fourth coil, a better defined spot is obtained. 


To transform this point to a coordinate system 
fixed to the head, with the origin halfway between 
left and right ears, the x-axis pointing to the nasion 
and y-axis pointing to the left ear and the z-axis is 
perpendicular to these axes, we use the following 
transformation: 
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Here y N , yL and y R are the positions of the nasion, 
the left and right pre aucular points, respectively. 




4. 


Results 


The localization of the coils with the method 
described above will be demonstrated with a few 
examples. The MEG system present at our lab is a 
whole helmet system of CTF Systems Inc. with 151 
gradiometers with a baseline of 5 cm. The noise 
level of the SQUIDS is about 10 fT/VHz. The MEG 
system is placed inside a three layered magnetically 
shielded room. In all experiments, the coils are 
energized during 2 s and the resulting magnetic 
fields are sampled at a rate of 1250 Hz. The 
energizing frequencies were chosen in such a way 
that their harmonics did not coincide with each 
other or with the 50 Hz power line frequency. 
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Figure 2. A comparison of the positions of the 
fourth coil relative to the first three fixed coils to 
the theoretical grid. Systematic shifts and rotations 
have been eliminated by fitting the theoretical grid 
onto the coil positions. 

4.1 General accuracy of the method 

The accuracy of the method was determined 
experimentally by fixing three coils on a piece of 
wood, on which a rectangular grid of 7 by 8 points 
with a 3 cm mesh size was depicted. A fourth coil 
with tip was used to indicate each grid point. The 
coils were energized with sinusoidal currents with 
frequencies of 5, 7, 12 and 17 Hz. The average 
distance between the coil set and MEG sensors was 
30 to 40 cm. 

For each grid point, all four coils were localized 
with the method explained in the methods section. 



To obtain a measure of the accuracy of the method, 
a theoretical grid of the same dimensions as the one 
designed on the piece of wood was fitted point to 
point to the positions of the fourth coil. Then the 
average distance between corresponding points of 
the theoretical and fitted grid was computed. 

Two experiments were performed, one in which the 
piece of wood was fixed and another in which it 
was moved several times in order to simulate 
patient head movement. In the first experiment, the 
average distance to the grid was .22 cm and in the 
second one, it was .24 cm. The distribution of the 
measurement points is shown in figure 2. No 
systematic increase of the errors near the grid 
boundaries can be observed. 



Figure 3. An example is shown of electrode 
positions, measured with the proposed MEG 
digitizer. Each electrode was measured four times 
after another, in order to quantify the 
reproducibility of the measurements. 

4.2 Electrode positions 

The method was also applied to determine 
electrode positions on a subject's head in a 
combined EEG/MEG experiment. In this 
experiment, MEG compatible electrodes were used 
mounted on a nylon EEG cap. Three coils were 
fixed on the nasion and the left and right pre 
auricular points. The fourth coil with tip was used 
to depict the electrodes. The subject was instructed 
to keep the head fixed during the 2 s needed to 
energize the coils and collect the data, but was 
allowed to move the head freely between 
measurements. 

The reproducibility was determined by measuring, 
with the same subject, all electrode positions four 
times after another. Of these electrode positions the 
average positions were determined (w.r.t. the 






Discussion 


nasion-ear coordinate system) and of each 
individual measurement the distance to the average 
position was determined. It was found that the 
average of these individual distances over all 
electrodes amounts to .23 cm. It was also found 
that the largest variability was at the mid-saggital 
plane. 

In figure 3, the resulting electrode positions are 
shown in the nasion-ear coordinate system, using a 
spherical projection. One observes that the cap was 
not placed exactly symmetrically on the head. For 
that reason, the symmetry line of the electrode grid 
is slightly rotated to the right. 

4.3 Headshape 

The MEG 3D localizer was also used to determine 
the shape of the head in the coordinate system of 
the MEG device. Figure 4 shows an example. The 
raw data set consisted of 248 points. Of this data 
set all points for which either nasion-left ear, 
nasion-right ear or left ear-right ear distances 
deviated more than .5 cm from the average were 
omitted, as well as all points for which the 
estimated coil callibration (=|m|) deviated more 
than 2 % from the average. After this omission 190 
high quality points remained. For visualization, the 
points were connected by Delaunay triangles and 
only triangles appearing on the foreground are 
plotted. The figure shows that the digitized points 
are nicely distributed over the head shape surface. 



Figure 4. The head shape measured with the four- 
coil method. 


5. 

Compared to the methods described in [1] and [3], 
our method has the advantage that simple coils can 
be used, instead of three orthogonal and equally 
calibrated coils having their dipole moments at the 
same position. One advantage is that construction 
errors, giving rise to systematic localization errors, 
are avoided. Another advantage is that not only the 
coil position is determined, but also its orientation, 
which makes is possible to attach a tip yielding a 
well defined localization point. 

The advantage of expressing the position of the 
fourth coil with respect to three other three coils 
which are localized each time the fourth coil points 
to a new spot is that the subject may move his head 
between the measurements. When a single 
Polhemus device is used instead, the subject’s head 
has to be fixed during the whole measurement 
session. 
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